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We theoretically study an exactly solvable Gamma matrix generalization of the Kitaev spin model 
on the ruby lattice, which is a honeycomb lattice with "expanded" vertices and links. We find this 
model displays an exceptionally rich phase diagram that includes: (i) gapless phases with stable spin 
fermi surfaces, (ii) gapless phases with low-energy Dirac cones and quadratic band touching points, 
and (iii) gapped phases with finite Chern numbers possessing the values ±4, ±3, ±2 and ±1. The 
model is then generalized to include Ising-like interactions that break the exact solvability of the 
model in a controlled manner. When these terms are dominant, they lead to a trivial Ising ordered 
phase which is shown to be adiabatically connected to a large coupling limit of the exactly solvable 
phase. In the limit when these interactions are weak, we treat them within mean-field theory and 
present the resulting phase diagrams. We discuss the nature of the transitions between various 
phases. Our results highlight the richness of possible ground states in closely related magnetic 
systems. 

PACS numbers: 71.10. Jm,71.10.Kt,75.50.Mm 



I. INTRODUCTION 

The study of zero-temperature properties of "frus- 
trated" magnetic systems has been reinvigorated in re- 
cent years with the development of powerful numeri- 
cal methodspH^ the discovery of new classes of exactly 
solvable models and the realization that intriguing 
topological ground states could occurP^l l n particular, 
there is a large class of exactly solvable quantum spin 
models known as Kitaev models^ which have revital- 
ized the study of exactly solvable spin-liquid systems. 
The appeal of this class of models lies in the relative 
ease in which Kitaev's original version can be gener- 
alized, spawning many variantsj^^l anc j its possession 
of many non-trivial properties^" even in the presence 
of disorder ESH2U Moreover, the exact solution of Kitaev 
models requires no more effort than solving the prob- 
lem of non-interacting particles moving in a static back- 
ground magnetic field. Nonetheless, the exact eigenstates 
of all Kitaev mode ls ar e non-trivial entangled many- 
body wavef unctions More specifically, their respec- 
tive ground states are examples of quantum spin-liquids, 
which are insulating quantum states of matter exhibiting 
no conventional long range magnetic order at zero tem- 
perature. Although no physical example of these models 
has been established in nature, there have been several 
experimental proposals for realizing themP^H Neverthe- 
less, the importance of this entire class of models lies in 
providing a strong case - by proof of principle - for the 
existence of exotic emergent phenomena in many-body 
quantum phases of matter such as quantum spin-liquids, 
as well as their utility as model systems for the study 
of non-trivial, non-perturbative emergent quantum phe- 
nomena in general. 

In this paper we study a new exactly solvable Kitaev 
model and examine its response to the inclusion of in- 
teractions that destroy the exact solvability. This will 



motivate us to analyze the order that is favored by these 
new interactions and their effects on the exactly solv- 
able ground state. The m odel is a Gamma matrix model 
(GMM) extensio n 1 11 ! 12 * 16 ! on the two dimensional ruby 
lattice. Previous studies on the ruby lattice have yielded 
a topological insulator," fractional quantum anomalous 
Hall states," and also topological anyons." Here we 
show that a GMM on the ruby lattice realizes a quantum 
spin liquid ground state with an unusually rich phase di- 
agram. The model is then generalized by the inclusion 
of Ising-like interactions that spoil the exact solvability 
of the model. We explore the interplay between these in- 
teractions and the many ground state spin liquid phases. 

The paper is organized as follows. In section II, we in- 
troduce the Hamiltonian and discuss its solution by map- 
ping onto a system of non-interacting Majoranas moving 
in an emergent Z 2 gauge field. In section III, we discuss 
the various ground state phases that the model realizes 
and present a phase diagram. In section IV we exam- 
ine the Ising order in one of the gamma matrix opera- 
tors which may be induced by including Ising-like inter- 
actions. In section V we treat these extra interactions 
within the mean field approximation and discuss their 
effects on the phase diagram. We then end with conclu- 
sions in section VI. 



II. GAMMA-MATRIX MODEL AND 
HAMILTONIAN 

We consider spin-3/2 moments located on the sites of 
a ruby lattice, shown in Figflja). Alternatively, one 
can think of each site as being occupied by a single 
spin- 1/2 which has an additional orbital degree of free- 
dom. At each site we introduce five 4x4 Hermitian 
Gamma matrices which satisfy a Cli fford al gebra locally 
{17, It} = 2S ab , with a,b = 1, . . . gSEHHl. These ma- 
trices also have a representation in terms of bilinears of 
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FIG. 1. (a) The ruby lattice. The ruby lattice can be viewed 
as an "expanded" honeycomb lattice with triangles replac- 
ing the vertices and squares replacing the bonds, (b) The 
couplings in our Hamiltonian Eq.|3| and ground state flux 
configuration of the Z2 gauge field. Note the J' couplings 
(not shown) are also defined on the same links as the J's. 
The arrows indicate the choice of gauge used for the My fields 
in {5). 



SU(2) spin-3/2 operators on each site i as follows, 
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r? = (s?) 
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The Gamma matrices may also be expressed as tensor 
products of Pauli spin-1/2 matrices {er^r"} on the site 



r 2 
r 4 
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? ®r?, (2) 

where cr? = = 1, the 2x2 identity matrix. In the spin- 
orbit interpretation ([2| , a acts on the real spin degree of 
freedom and r on the orbital degree of freedom. Thus, 
one can view the Gamma matrix model we study below 
as applying to both a spin-3/2 model and a two-orbital 
spin-1/2 model. In either case, there are 4 physical states 
per lattice site. 

The Hamiltonian we study is given by 
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TL — Ho + H5, 

h = Ja £ rjr 2 + J' A J2 r ' 5r 1 

(ij)eA (ij)eA 

+ Jo £ r; ! r) • £ rf rf , 

<«>eO («>eO 
H 5 =-J 5 J2 T l (3) 

where we define the operators Tf b — [Tf, T^]/(2i) in 
terms of the Tf given in either Eq.Q or Eq.fla). The 
couplings Ja , i an d </£> are taken to be positive 
and are link dependent as shown in Figfljb). For generic 
couplings, the Hamiltonian has translational and six-fold 



rotational lattice symmetry. With the interpretation in 
terms of spin-3/2 moments, it has global Ising spin sym- 
metry under 180° rotations about the z-axis, and time- 
reversal symmetry (TRS), although TRS will be sponta- 
neously broken in the ground state as we will show below. 

While the model |3| has unusual spin symmetries 
in terms of the underling spin and orbital degrees of 
freedom, its structure allows an exact solution. A 
key ingredient for the exact solvability of the model is 
the existence of conserved operators on every plaque- 
tte of the lattice which commute with the Hamiltonian.- 
These operators are the three (Wa), four (Wb), and six 
point (Wo) operators defined on the triangles, squares, 
and hexagons of the lattice (as their subscripts sug- 
gest). Explicitly, they are given by Wa = Tj 2 T^ 2 Tl 2 , 

Wo = rf rf rfrfr^, and Wa = Tf>T?T™r}*, 

where the sites are labeled counter-clockwise and for 
Wa the first link (ij) lies on a triangle. The model is 
then solved by introducing a Majorana representation of 
the T-matrices using six flavors of Majorana operators, 
U7-C"-i'; ! - fiS Ci, di}, at each site i: 

r? = Ti = iad h rf = itfd h (4) 

where the Majorana fermions have the important prop- 
erty that (£?)t = £?, c\ = a, and d\ = d t . The Hamilto- 
nian written in terms of these operators then reduces to 
that of two species of Majorana fermions, c and d, which 
are non-interacting and move in a background Z2 gauge 
field Uij, 

Zltij &iClj 

(ij)eA {ij)eA 
+ Jo £ iuijCiCj + Jq £ VUi%j did j 



{ij)eO 
Jr^y^jcidj. 



(5) 



tin are defined by My = —iQS 2 if r/ € A 



Here the fields u, tJ 

and = —i£f£j if ij G O. The Uij are in general quan 
turn fields with eigenvalues ±1 since ujj — 1, hence their 
identification with a Z 2 gauge theory. Moreover, we can 
simultaneously diagonalize H and } since they com- 
mute. We emphasize that the interacting Hamiltonian 
([3]) with the representation Q is reduced to the effec- 
tively non-interacting Hamiltonian ([5| because the Mis- 
behave as constants in each flux sector of the Wa, Wa, 
and Wq. 

However, the full Hilbert space spanned by the Majo- 
rana fermions is overcomplete, and a constraint must be 
enforced to ensure that the Clifford algebra of T opera- 
tors is satisfied. This constraint is expressed by the oper- 
ator equation A = -rjr?rf Ifl* = -i&&$&adi = 1 
and is enforced by the projector P — FJ ; [ 1+ 2 Pi ] , where 
the product is over all sites in the lattice. The original 
Hamiltonian is obtained through H = PHP. 

The Z2 gauge fields define fluxes <fi p via exp(itp p ) = 
Yljkep ^ u jk, where jk is taken counterclockwise on each 
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elementary plaquette p. These fluxes are then related 
to the conserved W's by W p oc Yiijep u iji where ij is 
also taken counterclockwise. Since the eigenstates of the 
Hamiltonian are also eigenstates of the fluxes, once the 
fluxes have been specified, the Uij — ±f's are uniquely 
specified up to Z2 gauge transformations, and the Hamil- 
tonian H is then diagonalized and projected to yield to 
eigenstates of %. The ground state is then determined by 
a many-body Majorana wavefunction minimizing the to- 
tal energy. Due to translational symmetry, % describes a 
bandstructure for the c and d fcrmions. It can be shown 
that the minimal energy configuration is the one where 
all the negative "eigenst ates" o f the effective band Hamil- 
tonian ([5]) are occupied! 5 1 11 1 12 1 

One finds that 4> p = ±ir/2 in the triangular plaque- 
ttes, and <p p = 0, n in the hexagonal and square plaque- 
ttes. Under time reversal symmetry, W p — > ±W P where 
— (+) is for triangle (hexagon and square) plaquettes; it 
follows that (j) p — ¥ —<j) p for triangle plaquettes while <f> p 
remains unchanged for hexagon and square plaquettes. 
Consequently, a ground state with a certain flux pattern 
{<fi P } spontaneously breaks time reversal symmetry. The 
ground state energy of a flux configuration must be de- 
generate with the flux pattern obtained from {<j> D } by 
changing (j) p — > — (f> p on all triangular plaquettes jEHU 

To determine the physical ground state, the flux con- 
figuration which minimizes the ground state energy needs 
to be determined. This was accomplished by first numer- 
ically determining the ground state energy for each flux 
configuration for the symmetric couplings Ja = J'^ = 
Jo = Jq = J and J5 = 0. The ground state flux 
with the least energy is depicted by FigQb). Gener- 
ally the ground state flux will be a function of J5, but 
for simplicity we will not consider additional flux con- 
figurations in this work because an additional term (de- 
pending on Wa, Wb, Wq) can always be included to fa- 
vor a particular configuration without destroying the ex- 
act solvability! 20 * 28 ^ By general arguments regarding the 
Majorana represe ntati on of the spin operators in the 
3/2 represent at ion P^H the spin-spin correlation functions 
of the ground state are identically zero beyond nearest 
neighbor sites. Hence, the ground state is also a quan- 
tum spin-liquid in these observables. 

With the gauge sector of T-L specified, we can solve our 
system by mapping the Majorana fermions to complex 
fermions a,i, defined by c; — ai + a\ and di — —i(a,i — a\). 
This puts our Hamiltonian into the form 



m 1 J 

-2J 5 J2(alai-^), (6) 

i 

where we have defined the quantities iufj — iUfjJA(O) 
and iufj — iv>ij J^r™ for ij € A(o). This type of 
Hamiltonian is solved by taking a Bogoliubov transfor- 



mation in momentum space. Taking the Fourier trans- 
form of (JTjj) and introducing the Nambu spinor as ^ (k) = 

(ai(k) a,2[k) ••• a>§{k) a\(— k) ■■■ aj(— k)j where 
aj(fc) is the Fourier mode of the fermion creation operator 
on the ith site in the unit cell, we obtain 

H= &(k)H(k)*(k), (7) 

k£C+ 

with the identifications 

H(k) _(h(k) A(k) \ 
H[k >- {A(k) -h T (-k)) ' 

hij{k)= i (u^(fe) + «£•(*)) +2J B 6ij, 
^ii{k)=^i(u^k)-u^k)), (8) 

where the w's are the Fourier modes of the corresponding 
it's defined above, and where C + is half of the Brillouin 
zone. We only sum over half of the Brillouin zone to avoid 
the double counting introduced by the Nambu spinor. 
Because H(k) = —H(—k), the eigenvalues of H(k) ap- 
pear in pairs {Ej(k), — Ej{— k)}, j = 1,2,..., 6, and our 
model is particle-hole symmetric. This is a consequence 
of our Hamiltonian being written in terms of the Nambu 
spinor, which satisfies ^](—k) = ^i(k); when we diago- 
nalize H(k) 1 half of the states are redundant. 



III. BAND STRUCTURE AND PHASE 
DIAGRAM 

For general couplings the ground state may be gapped 
or gapless. For the gapped phases, the band structure 
may possess a non-trivial (non-zero) Chern number. We 
will refer to these non-trivial gapped phases as Chern 
phases. The physical consequence of a non-zero Chern 
number is the appearance of chiral gapless edge modes 
in a system with a boundary. If the system has a fi- 
nite Chern number, it will exhibit a quantized thermal 
Hall conductance. We find these signatures of a non- 
trivial topological phase by calculating the Chern number 
numerically * 29 * 43 ! and diagonalizing a system with bound- 
aries (a "strip" geometry) to determine the existence of 
gapless edge modes. This is most conveniently done by 
taking a Fourier transform to moment um space along the 
"length" of the strip geometry! 5 1 11 * 12 1 

An example band structure is shown in Fig[2j which 
illustrates the dispersion of Majorana excitations E{k) 
as a function of the two planar momentum components 
k x ,k y . The state in Fig|2] is gapped about zero energy 
(below which all states are occupied and above which all 
states are empty) and has a non-trivial Chern number. 
There are six bands above and six bands below zero en- 
ergy, corresponding to the 6 sites in the unit cell of the 



FIG. 2. Bandstructure for the coupling parameters 
{J a, J'a, Jo, Jo. J 5}={1-0, 0.75, 1.0, 1.33,1.0} where qi = 
(2tt/3)(V3 - 1) and 172 = tt(1 - There are 12 bands 

which exhibit a redundant particle-hole symmetry and the 
energy spectrum is gapped with Chern number 2. 




FIG. 3. Spectrum of the strip geometry of the sys- 
tem as a function of the momentum k y , parallel to the 
edge. The system has Chern number -1 and couplings 
{J a, J'a, Jo, Jo,Js}={1-0, 0.6, 1.0, 1.67,0.2}. Two chiral 
edge modes that traverse the band gap at E — are clearly 
visible and each corresponds to a mode on a single edge. 



ruby lattice. Shown in Fig[3] is the spectrum of a strip 
geometry in the non-trivial phase with Chern number -1. 

When the ground state is gapless, it generically has a 
Fermi surface. Shown in Figs. [4] and [5] are band struc- 
tures with a Fermi surface realized in the ground state. 
In Fig|4] the Fermi surface is in the center of the first 
Brillouin zone (BZ), while in Fig|5]the Fermi surface is 
near the high symmetry K and if points of the hexago- 
nal BZ. Note that the Fermi surface is "connected" once 
the appropriate reciprocal lattice vectors are used to shift 
them together. At a phase boundary between two gapped 
phases with different Chern numbers, the ground state is 



-34 L * • * • • ,_ 

-qi q\ 

*x 

FIG. 4. (Color online) Bandstructure and shape of 
the Fermi surface (blue) at { J a, J'a, Jo, J'q, Js} = 
{1.0,1.0,1.0,1.0,1.4} and where gi = (2tt/3)(\/3 - 1) and 
q2 = 7r(l — l/v3). The Fermi surface pocket is located around 
Brillouin zone center. The dots are indicative of the zone 
boundary. 



gapless with Dirac nodes at a discrete set of points in 
the BZ. An example of such a band structure is shown 
in Figj6j which has four nodes in the first BZ (one at the 
zone center and three others at the inversion symmetric 
M points). 

In general, the phase diagram is complicated due to 
the large number of tunable parameters, even when one 
insists on translational, six-fold rotation and inversion 
symmetry. We report on a phase diagram resulting from 
an interesting selection of parameters which reveals the 
general richness of the model. We explored the phase di- 
agram for the specific case where the couplings Ja = Jo 
are held equal to a fixed constant J, and the J' couplings 
are allowed to vary in such a way that their product re- 
mains a constant; explicitly Jq x J' a = J 2 . This phase 
diagram is shown in Fig(7] and is exceptionally rich with 
many gapped and gapless phases. 

The regions of the phase diagram are classified into 
whether or not they are gapped or gapless, and whether 
or not the gapless phases possess a non-trivial Chern 
number (Chern phase). If gapless, the phases are classi- 
fied according to whether the gap closes at discrete Fermi 
points in the BZ or along a line or Fermi surface (FS). 
Often the former contains Dirac nodes on the BZ bound- 
ary or the zone center where the Fermi point may also be 
a quadratic band touching point. The Fermi points only 
occur on critical lines separating two Chern phases or a 
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FIG. 5. (Color online) Bandstructure and shape of 
the Fermi surface (blue) at {J/s, J' A , Jq, Jq, J5} = 
{1.0,1.0,1.0,1.0,0.55} and where gi = (2tt/3)(V3- 1) and 
52 = 7r(l — l/v3). The Fermi surface pockets are located 
around the K and K' symmetry points of the Brillouin zone. 
The dots are indicative of the zone boundary. 



at constant J' A /J = 1 where = J' A = Jq = Jq = J. 
This line contains an gapless FS phase in a line between 
0-7 % J5 < 2.0, and there are four NLs which cross this 
FS line at its endpoints. Along this line there are no 
nearby Chern phases in the phase diagram. This FS 
phase is non-critical in regards to tuning J' A and Jq, 
appears to be particularly stable, and is the FS phase 
with the largest volume in the phase diagram. 

Finally, we note that in the phase diagram transitions 
between Chern phases and FS phases involve critical NLs. 
This transition will lead to Fermi surface pockets forming 
around the nodes located at the points of high symmetry. 
Figures [4] and [5] are examples of such pockets. There 
are also transitions between Chern phases and FS phases 
which do not involve a NL. In this case the Fermi surface 
becomes gapped by a pairing (Cooper) instability as the 
system is driven into the Chern phase. 

In the next section of this paper, we will study the 
effects of the inclusion of new interaction terms which 
spoil the exact solvability of the model. These new in- 
teractions will also motivate the study of the Ising order 
in one the gamma matrices, namely rf . 



IV. ISING ORDER OF T 5 

In this section we consider the following extension to 




FIG. 6. Bandstructure of a nodal phase with four 
nodes and at parameters {J&, J'&, Jq, Jq, J5} = 
{1.0,0.58,1.0,1.73,2.62}. The nodes are located at the 
zone center and the three M points of the Brillouin zone. 



Chern phase and a Fermi surface. We will refer to these 
critical phases as nodal lines (NL). 

The phase diagram in FigjT] contains NLs with up to 
3 Fermi points. Since the NLs are obtained by solving 
a secular equation of a Bloch Hamiltonian which is ana- 
lytic, the NLs trace out smooth trajectories in the phase 
diagram. Moreover, the complicated manner in which 
they cross yields the many Chern phases in FigjTJ Most 
interestingly, we see a 3 node NL which closes in on itself 
to form an ellipse. 

Another interesting region is the high symmetry line 



H — Hq + %5 + "Wising) 

n Ising = -\J2^ 5 r (9) 

(ij) 

We also limit ourselves to case where A > 0. For our 
analysis it is convenient to take the spin-orbit interpre- 
tation of the model where we have 4 local degrees of 
freedom with basis states {|tt) j Itl) , lit) > III)}- In this 
basis, r| is trivially diagonal with eigenvalue +1 if the 
1/2— spins are aligned and eigenvalue —1 if they are anti- 
aligned. Thus, we can regard T 5 as an Ising spin, albeit 
one that is doubly degenerate for each Ising polarization. 
Then "Hising, as its name suggests, describes ferromag- 
netic nearest neighbor Ising-Ising interactions in the case 
where A > 0. It is also trivially true that ["Hi S ing,rf] = 
for all sites i, which implies that both operators are si- 
multaneously diagonal in the local spin-orbit basis. Thus 
in the "Ho = limit, the model is trivially exactly solv- 
able with eigenstates given by tensor products of local r| 
eigenstates. A typical eigenstate would have the form, 

IV') = <8> (on ltt> + ft III)) <g> (7, III) + Si lit)) , 

iei je.J 

ki 2 +ifti 2 = i7,f+i<y 2 = i, (10) 

where the N sites of the lattice are partitioned into dis- 
joint sets I and J with Ni and Nj sites each respectively. 
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FIG. 7. (Color online) Phase diagram of the spin-liquid at couplings J a = Jo = J an d J a x J'q = J 2 - The different 
colored regions correspond to either a gapped phase with Chern number as given in the legend or if gapless correspond to 
a Fermi surface (FS) phase. Lines which are boundaries between gapped phases are gapless phases with nodes or Fermi 
points at specific points in the Brillouin zone. The single node phase is gapless at the zone center k = (0, 0) and can be 
either a Dirac cone or a quadratic band touching point; the double node phase is gapless at the K and K' points with 
k = (±7r(y / 3 — l)/3,±7r(l — l/v3j), and the triple node phase is gapless at the three inversion symmetric M points, k = 
(0, ±7r(l — l/v3)) and k = (±7r(v3 — l)/2, ±7r(l — l/y / 3)/2). These lines are labeled according to the number of such nodes. 



The energy of such a state is then given by 

(H 5 + H Tsing ) |0) = {(m - Nj)J 5 - X(N }ink - L dom .)} |0) , 

(11) 

where JV un k is the total number of nearest neighbor links 
and Ldom. ^ s tne total length of all the boundaries be- 
tween the different "Ising domains" . The arbitrariness in 
the parameters a,/3,7,<5 on every site leads to a macro- 
scopic degeneracy of the state \ip); this stems from the 
local degeneracy of Tf . Also, an exact eigenstate of Tf 
will exhibit no fluctuations in the Tf observables. Hence, 
the ground states of "Hismg + %5 are trivial Ising ferro- 
magnets with order in (Tf) = ±1. The sign of the order 
parameter will depend on the sign of J5 if it is non-zero. 
Otherwise, if J5 = 0, the ground state will spontaneously 
break the Ising symmetry in r|. 

However, when Hq =/= the situation is no longer trivial 
and no longer exactly solvable. In particular, [%o,rf] 7^ 
suggests that it would be difficult to diagonalize the 
Hamiltonian in the spin-orbit basis or any basis where 
is diagonal. Nevertheless, one still has [Wp,%i sm g] = 
for all plaquettes p. This implies that the flux invariants 
are still good quantum numbers for both the GMM and 
Rising, and it is sensible to solve for the eigenstates within 



a given flux sector {Wp}. Writing %i S m g after performing 
the Majorana transformation makes this explicit, 

Wising = -4A J2 "'"•"/•■ + 2X J2 °i°< ~ 2NX - ( 12 ) 

Hence, Wising does not couple to the Z 2 gauge field de- 
grees of freedom and the fluxes are conserved. If we now 
consider again the limit where Hq — 0, then T-L describes a 
trivial insulating Hamiltonian with no kinetic terms but 
with local a-number conservation. One can then write 
down the exact eigenstates in this representation which 
is just any occupation state of the complex fermion cij for 
all sites i where the flux configuration can remain arbi- 
trary. This gives another explanation of the macroscopic 
degeneracy of \ip), which in the gauge field picture is due 
to the many degenerate {W p } sectors. Equivalently, the 
flux invariant operators are insensitive to whether or not 
the local 1/2-spins are aligned or anti-aligned; rather, 
they operate on degrees of freedom that are orthogo- 
nal. We can now also give a physical interpretation to 
the complex fermion a. Namely, its occupation parity 
Pi = (2a\ai — 1) is the Tf observable whose eigenvalue 
determines whether or not the local 1/2-spin are aligned 
or anti-aligned. Also, the a occupation numbers do not 
fluctuate in an eigenstate, which is consistent with the 
non-fluctuating behavior of the Ising order parameter. 
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Now we return to considering general H. Note that 
Ho neither commutes with Tf nor its sum total X^^f. 
This leads to the non-conservation of the total a-particle 
number in the ground state, and we can associate 
this non-conservation with the pairing terms in the a- 
representation of the problem. Nevertheless, the product 
of parities as defined by V — ]X Pi = Yli rf is conserved 
since Ho,H^ and "Hising all individually commute with 
it. In fact V implements a global Z2 gauge transform, 
en <— — cii, a\ <— — <zj on every site i. This is merely 
another manifestation of the conservation of the total a- 
number modulo 2 by H. 

Next we will consider the opposite limit where "Hising + 
H$ = and argue that the exact ground states of Ho 
will not possess any Ising order in T 5 by themselves, but 
that the exact ground states of the more general GMM 
with J5 7^ will. Moreover, in the limit of small J5, 
these ground states are adiabatically connected. Given a 
fixed flux sector {W^}, and given the existence of a non- 
degenerate ground state of Ho which we denote by |fi), 
we will show that (fi|rf \Q) = for all i. 

First consider equation ^ in the limit J5 = where 
the gauge fields Uy have been fixed. Define a linear uni- 
tary particle-hole conjugation operator C by its by con- 
jugation on a, and a\ and its action on the a-number 
vacuum |0) as follows: 



C|0) 



Ca[C^ = a,i, 



|0). 



e i6 a] . 



(13) 



9 is a phase which we can fix by demanding that C 2 |0) = 
|0). For TV even this leads to 9 — ^f-. These relations 
then totally specify C and also imply the relations = C 
and C 2 = 1. Note that [Hq,C] = 0, implying that if is 
non-degenerate, then C|0) can only differ from |f2) by a 
phase. Thus, (fi|ato*|fi) = aia\C\tl) = (fl\ ai a\\n). 



a]a\ com- 
1 j 



Using {a^aj} = 1, we then conclude that (f2|ajcij|f2) = 
1/2 and (f2|rf |fl) = 0. Since a\a,i is Z 2 gauge invariant, 
we expect this to hold true even after projection. One 
could also heuristically argue that because J5 couples to 
a\a.i like a chemical potential, the ground state will have 
a particle-hole symmetry and is the half-filled Fermi sea. 
However, the presence of pairing terms ataj 
plicates this argument. 

For the model that was solved in Sec. II where J5 7^ 
generally, (rf } may take non-zero values in (0, 1] when 
in the ground state. In this more general situation, the 
presence of breaks the symmetry under C conjuga- 
tion. In addition, since the ground states of the GMM 
were shown not to undergo a phase transition as J5 is 
tuned from zero, we conclude that the exactly solvable 
ground states of the GMM may exhibit Ising order in the 
rf spins. That is, the exactly solvable phase and Ising or- 
der are not mutually exclusive. Hence the ground state is 
not a spin-liquid in these observables. But fluctuations 
in rf remain non-trivial in the exactly solvable phase. 
This poses the interesting question as to what extent an 




Chern Number: 



-3 
_2 

-I 

+ ] 
FS 



FIG. 8. (Color online) Phase diagram with the Ising per- 
turbation treated within the mean field approximation. The 
couplings are { Ja, J'a, Jq, Jq} = {1-67, 1.0, 0.3, 1.0}. 



exactly solvable GMM ground state is adiabatically con- 
nected to the trivial Ising states such as the one described 
by above. 

In the limit where Ho is treated as a perturbation to 
H5 + Rising, the weakly fluctuating Ising order r| will 
strongly renormalize the J5 coupling as seen by the a 
fermions, which can be seen from equation ( 12 ). Thus, we 



can expect this limit to be equivalent to a large J5 limit 
and to be described by a low energy effective Hamiltonian 
derived from perturbation theory. Schematicall y, the 
Hamiltonian can written in terms of W p operators] 5 * 11 * 44 ! 



Hn 



Yl a P w p 



(14) 



where {a p } are coupling constants and the dots represent 
higher order terms. Such an effective theory will break 
the degeneracy among the different flux sectors and favor 
certain flux sectors over others. However, time reversal 
symmetry is still preserved and must be spontaneously 
broken by the ground state. If J5 = 0, then the ground 
must also spontaneously break the Ising symmetry of T 5 . 

In the opposite limit where Hi s i ng is the perturbation 
to Ho + H§, the situation is more complicated as the 
ground state of the exactly solvable phase is in general 
non-trivial. We try to address this question in the next 
section of the paper where Hi slng is treated as a perturba- 
tion to the exactly solvable GMM within the mean-field 
approximation. 



V. MEAN FIELD APPROXIMATION OF THE 
ISING PERTURBATION 

In this last section we study the effects of Hi sm g as 
a perturbation to the exactly solvable Ho + H§ within 
the mean-field approximation. Performing a mean-field 
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FIG. 9. (Color online) Phase diagram with the Ising pertur- 
bation treated within the mean field approximation. The cou- 
plings are { J A , J' A , Jo, J'q} = {1.0,1.0,0.2,0.2}. The phase 
boundary separating the two trivial phases is a first order 
transition line where (rf } changes discontinuously. 

decomposition with respect to the order parameter (r|) 
results in the following mean-field Hamiltonian, 

n Ising « H MF = -4A^(r 4 5 )( a |a 4 -i/2)+A^(rf)(r 3 5 ), 

(ij) (ij) 

(15) 

where the mean field (rf) has to be determined self- 
consistently. Motivated by the fact that the Ising cou- 
plings are ferromagnetic, we make the ansatz that (rf ) is 
uniform across the entire lattice. At the mean-field level, 
the effect of the interaction "Hi s i ng can already be seen to 
renormalize J 5 with jf = J 5 + 2\(T 5 }. 

Within this approximation, we determined the phase 
diagrams for varying J5 and A with fixed J and J' cou- 
plings and fixed flux sector. Two phase diagrams are 
shown in FigjS] and Fig(9] We note that generically the 
topological phases are stable against weak "Hising pertur- 
bations. At larger A, however, the energetics always pre- 
fer the trivial gapped phase which has zero Chern number 
and large (r 5 ). Hence, large A eventually renormalizes 
J5 to larger effective values. This is consistent with the 
large J5 and large A trivial phases being adiabatically 
connected to each other. 

Focusing along the line at constant A = and increas- 
ing J5, we see a series of gapped and gap less phases 
which eventually ends with a trivial phase at large J5. 
This is expected from the results of the previous section, 
but since J5 acts like an external field, the Ising order is 
(rf) ^ on this line. When A is tuned to greater val- 
ues, the volume of these intermediate phases diminishes 
continuously, eventually giving way to the trivial phase. 
This confirms that the trivial Ising phase is adiabatically 
connected to the trivial large J5 phase. 

In FigjHJ the gapped phase at J5 = appears to be the 
most stable with respect to increasing A. However, for 
the couplings shown in Figj9j the J 5 = phase gives way 
to a trivial gapped phase (Chern number zero) which is 



not adiabatically connected to the trivial phase at large 
J5 . This phase undergoes a first order phase transition to 
the larger J 5 trivial phase as shown in Figj9j Thus, the 
Ising interaction may choose to stabilize certain gapped 
phases over the default J5 = A = phase. The stabilized 
phases are highly dependent on the particular values of 
the other coupling constants. 

At constant J 5 > 0, increasing A may stabilize certain 
phases which would typically require larger J5 values. 
For example, in FigjH] along the line J5 = 0.75, the sys- 
tem is driven to a Fermi-surface phase with increasing A. 
Other examples of such transitions driven by the Ising 
interaction may also been seen in Figs. [8] and [5J Hence, 
the Ising perturbation treated at mean-field may stabilize 
certain exotic phases which are otherwise only accessible 
by intermediate values J5. 

In summary, these results show that with increasing 
A, the ground state is eventually adiabatically connected 
to the trivial Ising ordered state of Sec. IV. However, in 
general, various phase transitions involving other exotic 
phases must occur before this happens. Moreover, this 
depends crucially on the initial (A = 0) exactly solvable 
phase. 



VI. CONCLUSIONS 

In this work we studied an exactly solvable T-matrix 
generalization of Kitaev's original spin-1/2 model on the 
ruby lattice shown in Fig[T] Our model can be inter- 
preted as describing either a spin-3/2 system or a double- 
orbital spin-1/2 model. We have shown that the ruby lat- 
tice, with its many possible tunable parameters, exhibits 
a rich phase diagram with many exotic phases realized 
by same Hamiltonian. We have found gapless phases 
with Fermi surfaces and Fermi points, as well as gapped 
topologically non-trivial phases. Next, we studied the 
effects of including Ising-like interaction terms that de- 
stroy the exact solvability of the model. We analyzed 
the ordered phase that these terms lead to and argued 
that the ground state in the large J5 limit of the exactly 
solvable model is adiabatically connected to it. We then 
analyzed in detail the general model near the limit of 
an exactly solvable phase and treat the Ising terms as 
perturbations within the mean-field approximation. We 
derived phase diagrams which show that the Ising in- 
teractions will eventually favor the trivial Ising ordered 
phase of the model over the other more exotic phases. We 
also confirm that the large J5 exactly solvable phase is 
adiabatically connected to the trivial Ising ordered phase. 

While this study has focused on a particular exactly 
solvable spin model and special perturbations on it, con- 
nections previously drawn to other topological phases in- 
dicate that its results are rather broadly applicable.^ 
Moreover, a number of future directions for further study 
suggest themselves: (i) Considering anti-ferromagnetic 
Ising interactions (A < 0) to determine if the interplay 
between the exactly solvable GMM and the Ising in- 



9 



teractions might be different, (ii) The effect of disor- 
der in interacting many-body quantum systems remains 
a key open problem. The model we study here with 
its unusually rich phase diagram provides an opportu- 
nity to study disorder effects within a well-controlled 
models (iii) Doping magnetic systems is one route to 
high-temperature superconductivity. To date, no work 
has been reported on the doping of Gamma matrix gen- 
eralizations of the Kitaev model, and only a few stud- 
ies exist on the original model The richness of the 
model we study here could shed light on doping mag- 



netic systems in a much more general context and could 
ultimately help guide the discovery of strongly correlated 
materials exhibiting high temperature or unconventional 
superconductivity. 
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